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Abstract. We show that the refined Donaldson-Thomas invariants of C 3 , 
suitably normalized, have a Gaussian distribution as limit law. Combinatori- 
ally these numbers are given by weighted counts of 3D partitions. Our tech- 
nique is to use the Hardy— Littlewood circle method to analyze the bivariate 
asymptotics of a q-deformation of MacMahon's function. The proof is based 
on that of E.M. Wright who explored the single variable case. 



1. Introduction. 

In [7] physicists suggested that the (IIA) string theory associated to a Calabi-Yau 
threefold X should produce an algebra of BPS states: 

Hbps = 

7er 

where here the (charge) lattice T is identified with the even cohomology of X, 
r = (B^ =0 H 2l (X, Z). Moreover each individual vector space % 1 should have an 
additional Z-grading coming from a symmetry in the little group Spin(3) [BJ. 

Mathematically we consider the cohomological Hall algebra [H] as giving the al- 
gebra of BPS states on X = C 3 . In this case T = H e (X,Z) = Z and the 7 = nth 
piece is given by the critical cohomology of Hilb"(C 3 |^] Moreover each of these 
vector spaces has a cohomological Z-grading. The Betti numbers of these graded 
pieces are know as refined DT invariants. These numbers are dependent on the 
singularities of Hilb"(C 3 ). For example they do not satisfy a Poincare duality. 

However, on a recent visit to EPFL, T. Hausel shared with me the output of a 
computer experiment. He conjectured that the refined DT invariants, suitably 
normalized, would have a Gaussian distribution as limit law, i.e. for large n plot- 
ting the Betti numbers against cohomological degree should give the bell curve of 
a Gaussian distribution (cf. [18J). The goal of this paper is to prove that conjecture. 

In fact this proposal is entirely combinatorial. The Hilbert-Poincare series for the 
cohomological Hall algebra, computed in [2J, equals M 3 (t, q 1 ^ 2 ) where t gives the 
T-grading, q gives the cohomological grading, and 

TCI— 1 _^ 

M s (t,q) = n n 1 _ a 8+2k+i- mtm - 

m>l fe=0 ^ 



In addition we add a single D6-brane filing the C 3 or mathematically a framing. 
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Expanding this series gives an explicit formula [14] for the t n coefficent 




7rhn 



where the sum is over all plane partitions of n which we now explain. 

A plane partition is given by a two dimensional array of positive integers in 
the first quadrant of 1? that are weakly decreasing in both the x, y directions. In 
this way plane partitions are a generalization of ordinary (line) partitions jlj. The 
analogue of the Young diagram in this situation is a stack of three dimensional 
boxes in Z> . Such a collection gives a plane partition if and only if the stack is 
stable under the pull of gravity along the (1, 1, 1) axis. For example let tt be the 
plane partition given alternatively as 





2 2 10 

3 3 2 

4 3 2 1 

5 3 2 1 1 




Figure 1. integers Figure 2. boxes 

Here the total sum of the integers/boxes is \w\ — 35 and we say that the partition 
has size 35. The statistics appearing in the above formula for refined DT invariants 
are, 

w + w = ^2 ni >3 > w - w = X! ^ ' w ° w = z2 

i<j i>j i 

Considering the set V n of plane partitions of size n as a sample space with 
uniform measure we define three random variables 

X£ > X n , X n :V n ^M. 

given by w+/n 2 ' 3 ,W- /n 2 ' 3 , and wq/u 2 ^ 3 . Our main result is the following: 

Theorem 1.1. The distribution of the random vairable 

5 ■ X°+ X+ - X n 

for large n has the Gaussian distribution as limit law with 

ft = 5C(2)/(2C(3)) 2 / 3 and a 2 = 1/(2C(3)) 1 / 3 . 

2. Setup. 

First we split the problem into two parts, one of which has already been solved. 
Straight away we see that the covariance of X^ and X+ — X~ is zero due to 
symmetry 

E((X°- Mx o)(X+-X-- Mx +_ x -)) = E((X° - » x o)(X+ - X-)) 

= E(X° n (X+-X-))-E( Px o(X+-X-)) 
= E(X°X+ - X° n X-) - Mx o (E(X+ - X-)) 
= 0-0. 
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The following is a result of E.P. Kamenov and L.R. Mutafchiev: 

Theorem 2.1 ([IT]). Let a = C(2)/(2C(3)) 2 / 3 and b = ^/l/3/(2C(3)) 1/3 . Then as 
n — > oo we have 

X°~Af(a,bn~^ In 5 n). 

Now the sum of two Gaussian random variables is again Gaussian with mean 
the sum of the means and variance the sum of the variances plus a covariance term. 
Since the above covariance was zero, Theorem |1.1| will now follow from the result 
just mentioned together with: 

Theorem 2.2. Let c = 1/(2C(3)) 1/3 . Then as n -> oo we have 

X+ -X- ~JV(0,c). 

To prove this result we use the method of moments. That is we show that the 
limiting distribution has the same moments as a Gaussian random variable with 
variance 1/(2C(3)) 1 / 3 . Specifically we will show that in the limit 



¥ (( Y + _ Y -\ k \ - J if fc is odd, 

n > ) \ (fc-l)!!(2C(3))- fc / 6 if fc is even. 



Consider the generating series 

m— 1 1 

m>0 fe=0 q tt 

and let p n (q) be the coefficient of t n then we have 

where d = q-^. Notice by symmetry this already implies that all the odd moments 
vanish. The method of proof given in the next section follows the proof of E.Wright 
[20j who provided an asymptotic formulgj^j for the number of plane partitions of n 

PnW\ q =l 2ll/36 /3^ n25 /36 e 



Wright's proof in turn generalized the pioneering work of Hardy and Ramanujan [5] 
who first applied the Hardy-Littlewood circle method to get an asymptotic formula 
for the number of ordinary partitions. Using this method in the next section we 
will show that 



when fc is even. This gives the correct moments and shows that X+ — X n has a 
Gaussian limit law as promised. 



2 Typo warning: Wrights formula on page 179 is missing factor of \/3 found at the end of his 
proof on page 189. 
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3. Proof. 

As explained in the previous section we are going to use the Hardy-Littlewood 
circle method to estimate the coefficients in the generating series 

M k (t) :=d k M {t,q)\ q=1 . 

Given any function A(t) — X) n >o a «^™ analytic on the interior of the unit disk we 
can compute the nth coefficient in its MacLaurin series using the Cauchy formula 

1 



2ni 



r n ^ 1 A(t)d* 

Cm 



where Cn is the circle or radius e" 1 ^ . The idea of the circle method is that by 
understanding the singularities of A on the unit circle one can approximate this 
integral by an integral over a small subarc of the circle when N,n^$> 1. 

For example when A(t) = M(t) = M(t, q)\ q —i is MacMahon's function then letting 
t = e z Wright defines the major arc C' N to be the points such that im(z) < 1/N 
and the minor arc C" N to be the remaining points on the circle Cn- In some sense 
MacMahon's function is most singular at z — and so the integral over Cn is well 
approximated by that over the small arc C' N . The following two Lemma^j make 
this precise and will be very useful to us later: 

Lemma 3.1 (E.M.Wright [5D] Lemma I). There exists constants N ,K so that for 
all N > No and t = e z G C' N along the major arc we have 



M(t)-e- 4 lV' u e 



C'(l) ,1/12^(3)/^ 



< A'7V- 1 /12^2 e C(3) 



AT 



Lemma 3.2 (E.M.Wright 20J Lemma II). Given any e > there exists an No 
such that for all N > Nq and t 6 C N along the minor arc we have 

\M(t)\ <e (C(3)-l/2+ 6 )JV 2 . 



Using Lemma 3.2 one shows that the integral along the minor arc is relatively 
small. Then using Lemma |3.1| the integral along the major arc can be approxi- 
mated using the curve of steepest decent. This gives Wright's asymptotic formula 
mentioned earlier [2D] and illustrates the idea of the circle method. 

Recall, we are specifically interested in computing the coefficients of the series 
Mk(t) — d k Mo(t,q)\ q= i defined at the start of the section as a means to com- 
pute the moments of the random variables X£ — X~. Let us write this series as 
M k {t) — F k (t) ■ M(t) then by Wright's two lemmas above we have a good under- 
standing of the singularities of the factor M{t) it remains to analyze F k (t). 

Example 3.3. Computing F 2 (t). Let us differentiate Mo(t,q) twice using d = 
this gives 



im>l fc=0 



M (t, q) 



■^Wright's first Lemma is more refined than this. We mearly extract his leading order approx- 
imation suitable for our purposes. 
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2 

■M (t,q). 

,m>l k=0 

Then setting q = 1 we get 

m— 1 

E E I • M ® = r I E I ■ w 

i>l fc=0 

deducing that 




(1 - i m ) 2 ' 
m>l v ; 

In the next two Lemmas we analyze the behavior of Fk{t) along the major and 
minor arcs. 

Lemma 3.4. Given k even, there exist constants No, K such that for all N > Nq 
and t — e z G C' N along the major arc we have 



F fc (i)-(fc-l)!!(2C(3)) fe / 2 z- 2fc 



< KN 2k ~ 2 . 



Proof. First let us consider the case k = 2 used to compute the variance. By what 
we saw in Example |3.3| above 

m(m 2 — l)e mz 



F, 



w 3 4^ (1 



(1 - e mz ) 2 



By using the Mellin transform e T = ^ j°_ %0 ° T(s)r s ds we can express this 
function as an integral, 

n 1 v-^ m(m 2 - l)e mz 1 \ - \ - . 2 

m>l V ' m>l i>l 



6 i/ . rwEE ,,,, 



m(m — l)i 

as 



1 

6z7T 



CT + ZOO 

r( 5 )z- fl c(s - i)(C(«- 3) - C( s - i))ds. 



(7 — TOO 



This integrand has a double pole at s — 2 and a simple pole at s = 4 coming from 
the Riemann zeta function while the gamma function contributes simple poles at 
s = 0, —1, —2, .... Doing the residue calculus we see that 

where 7 is Euler's constant coming from the Laurent expansion of the zeta function 
about 8 = 1. 

Computing the higher order derivatives is essentially an application of Wick's 
theorem. Indeed if k — 2r then we have 

F k=2 r(t) = (k-l)ll(F 2 (t)Y + G k (t) 

where the combinatorial coefficient (k — 1)!! comes from all possible ways of pairing 
the 2r differentials d 2r . Using the product rule for differentiation we see the first 
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term appearing. The remaining terms in G k (t) come from the other terms gener- 
ated in the product rule. Computing their Mellin transform shows the poles they 
contribute are of order at least two less. □ 



Lemma 3.5. Given k even, there exists positive constants No, C k , A k such that for 
all N > No and t = e z G C'^ along the minor arc we have 

\F k (t)\ < C k N A K 

Proof. Looking at the definition of F k (t) we see that ultimately it can be written 
as a finite sums and products of series like 



4-? a 



n a t r ' 



n>l 



(1 -t n ) b 



where a, b, c € N. Each of these sums can be bounded like 



n a t Tl 



E ll L. 
(1 - t n f 



< 



n a \t\ 



n>l v 1 1 ' 



using the Mellin transform as in the previous lemma gives bounds on this sum like 
C a ,bN Aa < b where \t\ = e~ x / N and A a t,,C a .b are constants depending only on a,b. 
In total this gives the polynomial bounds claimed. □ 

Now we have got bounds on the series M k {t) along the major and minor arcs we 
are ready to estimate the Cauchy integral. Let us define the two quantities we are 
interested in computing 



ji = 

k,n 



J" 
1 h 



1 

2ni 

1 
2iri 



r 



( AT 



C' 



L M k (t)dt, 



v M k {t)dt. 



From now onwards we choose N = (n/2((3)) 1 ^ 3 so that n — 2£(3)iV 3 . Now we 
can get a bound on the integral along the minor arc that is sufficient to show this 
integral is negligible: 

Lemma 3.6. Given k even then for all e > there exists positive constants No, K, A 
such that for all N > Nq we have 



Proof. Using Lemmas |3.2| and |3.5 | we get 

F k (t) ■ M(t) -t-^dt 



< 



1 

2~TTi 

|^ • sup(|F fc (t)|) • sup(|M(i)|) • supd*-"" 1 !) 
Zm\ nn n n rn 

I 1>jv °JV °JV 

= (C(3)"l/2+£)W 2 . J-n-^N- 1 



substituting n = 2((3)N 3 gives the result. 



□ 
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From now all that remains is to estimate the integral I n k . Combining Lemma 
3. II and Lemma 13^41 we have 



In,k = ^ [ F^M^r^dt 
2m J c , n 

= e f(-P (fe-i)!! rW /12 cv +2mN3 
2m (2C(3)) fc /6 J (1 _ i)/N 

+ 0(JV- 1+2fc - 1 /12-2 e 3C(3)iV^ 

" (2C(3))*/e 2izi 

+ 0(7v- 1 + 2 ' i; - 1 / 12 - 2 e 3 ^ 3 ) Ar2 ) 

where we set v = Nz. 

Notice the prefactor in the above expression is essentially the term we are looking 
for. From now on we work with the integral 

p' k := — f 1+l v -^+i/i2 e am 2 (2v + v^) dv 

basically it will be enough to show that, for large N, this is independent of k. Then 
in the limit we have I' n k = ^c^))^/ 6 ' ^«,o- We are able to achieve this by localizing 
the integral P' n k to an even smaller arc using the method of steepest descents. 

In the exponent of the integrand we have the function 2v+v~ 2 . Roughly speaking 
the integrand will be largest when this exponent is real. The curve of steepest 
descent is defined to be the real curve given by im(2v + v~ 2 ) = 0, specifically 
taking v = X + iY we have 

(X 2 + Y 2 ) 2 =X. 

This is the closed curve C meeting the lines X = at (0,0), X = 1 at (1,0), 
X = -Y at D = (2-1,-2-1), and X = Y at E = (2"l, 2~§) as seen in Figure |] 
As in [20] we consider an alternative integral along this curve C rather than along 
the straight line FG. 

Making a branch cut from to — oo along the real axis we consider the value of 
v 1 / 12 which is real and positive at v = 1 and take the contour for C parameterized 
in the anti-clockwise direction. Following Wright we define 

= e C(3)JV2(2 " + "' 2) ^d J' Kn = J c Uv)dv. 

On the straight lines EG, DF and along the arcs OE and OD we have good 
bounds on £ k (v). Since here Re(y- 2 ) = (X 2 — Y 2 )/(X 2 + Y 2 ) 2 < setting Y — rX 
gives 



v -2k+l/12 e C(3)N 2 (2v+v- 2 ) 



< X- 2k (l + r 2 )- k e Ci3)N yt^ X ~ 2 )e 2C( - 3)N2x . 



When r > 1 this tends to zero as X — > 0, and when r = 1 along the lines EG, DF 
there are easy bounds. In summary, we have a bound |f fe (t)| < Ke 2 ^ N2x along 
these contours OE, OD, EC, and DF. Using contour integration to compare the 
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Figure 3. Curve of steepest descent. 



original integral P' n k to the new integral J' n k along the curve C we see that 



















\pt _ V 1 < 
\ n,k ° n.k \ ^* 


/ £,k{v)dv 

J F 


+ 


f 

Jd 


+ 


f 

Jo 


+ 


f 

Je 



this allows us to integrate along the curve of steepest descent instead. To parame- 
terize this curve we choose 



t 



v- 1 



(2« + l)a 



so that t 2 — 3 — 2v — v 
line 



Now the problem transforms to an integral over the real 



J'u = e 3C(3)Af2 / Xk {t)e~^ )N2t2 dt, with Xk (t) 



v -2fc+l/12 dv 



2-ni 



dt 



The most serious piece of this integral is located about t = i.e. v = 1. To under- 
stand the behavior here we take a convergent power series Xk(t) — J2m=o a mt m in 
a small neighborhood of t = 0. Next observe that for some constant K we have 



dv iv 2 (2v + 1) 2 
dt ~ 1 + v + v 2 



w -2*+X/i2 dv 
2ni dt 



,-2fc| 



< K\v~ 

on all of the real line. Moreover for some possibly larger K we have 



\v- 2k \ = \2,-2v~t 2 \ k <Kt 2k 

on the compliment of the above radius of convergence about t — 0. All in all we 
get that 

2fe+l 

Xk (t)- J2 a ™ <m Kl 

m=0 



2k+2 
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over the whole real line. Finally our integral is approximated by 

2fe+l 



where 



/oo 
t m e -am^ dt + Mk 
., ., -oo 



\M k \ < K / t 2k+2 e-^ N2t2 dt < 



jV 2fe + 3 ' 

By symmetry all the above odd integrals are zero. The even ones are given by 

t 2m e C(3)N 2 t 2 dt = F ( TO + 2) 

(C(3)7V 2 ) m +5' 

So to get the leading order asymptotics we need only the constant term ao = 
in the expansion of \k (t) . In particular to leading order there is no dependence on 
k as we claimed earlier. Substituting N = (n/2£(S)) 1 ' 3 gives the formula described 
at the end of Section 2. 



4. Final Remarks. 

Some remarks about the asymptotics of DT invariants coming from this investi- 
gation and math/physics literature. 

Asymmetry. Since the refined DT invariants are given by M 3 (t, q 1 ^ 2 ) their distri- 
bution is shifted by the trace of the plane partition coming from 3X® . This shifting 
is relevant to the refined topological vertex in physics flu] . This was discussed in 
[13j where a minor discrepancy was noticed with calculations in [5] for the case of 
the refined DT invariants of the resolved conifold singularity. 

Dimension of moduli space. In his lecture notes on quiver moduli M. Reineke 
describes a conjecture on M. Douglas on the asymptotic growth of Euler numbers 
of spaces of representations of Kronecker quivers |17j . Reineke proposes a general- 
ization of the conjecture would give 

ln(x(M d (Q))) ~ C Qy /dim(M d (Q)) 

where Md{Q)) is a suitable smooth model for the quiver moduli, d is a large di- 
mension vector, and Cq is an interesting constant to be determined. 

For sheaves on a Calabi-Yau threefold the moduli spaces will be singular and 
not of the expected dimension. However in our case we do know that 

Cc 3 Vdim(Hilb n (C 3 )) < In (|# vir (Hilb"(C 3 ))|) < C+ v / dim(Hilb"(C 3 )) 

for large n, where # vir is the virtual Euler number or numerical DT count for this 
moduli space, and C^ 3 ,C^ 3 are constants. A proof of this follows from Wrights 
theorem |20j for the virtual Euler number and Briancon and Iarrobino's asymp- 
totics for the dimension of Hilb™(C 3 ) Geometrically this relationship between 
the virtual Euler number and dimension of the moduli space seems like a strange 
coincidence specific to C 3 ? 
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Orbifold: C x [C 2 /Z2]. In the case of numerical DT invariants, Panario, Rich- 
mond, and Young investigated the orbifold C x [C 2 /Z 2 ], here the charge lattice is 
two dimensional and one studies the bivariate asymptotics of colored partitions see 

US]. 

BPS black holes. A major goal of string theory is to unify quantum mechanics 
with Einstein's general relativity. In the 90's Strominger and Vafa showed that 
indeed the topological string described the physics of some BPS black holes in a 
certain limit of the string coupling constant |19j . Later with Ooguri they conjec- 
tured that in a certain limit the entropy of BPS black holes should be determined 
from the square of the topological string partition function [15]. 

Recently in the case of D6 - D2/D0 states this conjecture has been developed 
by Denef and Moore [I]. If we let X be a Calabi-Yau threefold and I(X,(3,n) be 
the moduli space of ideal sheaves with Chern character (1, 0, — n), their paper 
speculates that 

Um An (ln(|# v "(/(Y,A 2 /3,A 3 n))|)) 
7->oo I hi(A) 

By Wright's theorem this is true when (3 = 0. In [8 physicists checked this con- 
jecture for the quintic threefold. However as it is a hard problem to compute the 
higher genus Gromov-Witten invariants in this case these checks are only partial. 
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